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ABSTRACT 

We construct merger trees from the largest database of dark matter haloes to date 
provided by the Millennium simulation to quantify the merger rates of haloes over a 
broad range of descendant halo mass (10 12 < M < 10 15 M Q ), progenitor mass ratio 
(10 -3 £ ^ l)i an( i redshift (0 ^ z < 6). We find the mean merger rate per halo, B/n, 
to have very simple dependence on Mq, £, and z, and propose a universal fitting form 
for B/n that is accurate to 10-20%. Overall, B/n depends very weakly on the halo 
mass (tx A/q' 08 ) and scales as a power law in the progenitor mass ratio (oc £~ 2 ) for 
minor mergers (£ < 0.1) with a mild upturn for major mergers. As a function of time, 
we find the merger rate per Gyr to evolve roughly as (1+z)"™ with n m — 2 — 2.3, while 
the rate per unit redshift is nearly independent of z. Several tests are performed to 
assess how our merger rates are affected by, e.g. the time interval between Millennium 
outputs, binary vs multiple progenitor mergers, and mass conservation and diffuse 
accretion during mergers. In particular, we find halo fragmentations to be a general 
issue in merger tree construction from iV-body simulations and compare two methods 
for handling these events. We compare our results with predictions of two analytical 
models for halo mergers based on the Extended Press-Schechter (EPS) model and the 
coagulation theory. We find that the EPS model overpredicts the major merger rates 
and underpredicts the minor merger rates by up to a factor of a few. 



1 INTRODUCTION 



In hierarchical cosmological models such as ACDM, galax- 
ies' host dark matter haloes grow in mass and size primarily 
through mergers with other haloes. As the haloes merge, 
their more centrally concentrated baryonic components sink 
through dynamical friction and merge subsequently. The 
growth of stellar masses depends on both the amount of 
mass brought in by mergers and the star formation rates. 
Having an accurate description of the mergers of dark mat- 
ter haloes is therefore a key first step in quantifying the 
mergers of galaxies and in understanding galaxy formation 
and growth. 

Earlier theoretical studies of galaxy formation typi- 
cally relied on merger trees generated from Monte Carlo 
realisations of the merger rates gi ven by the analyti- 
cal extended Press- Schechter (E PS; lLacev fc Cold " 
Bond et al.l 119911 ) m odel ( e.g. iKauffmann et alj 



1993; 



1993; 



Somerville fc Primackl 1 19991 ; ICole et all l2000t ). Some re- 
cent studies have chosen to bypass the uncertainties and 
inconsistencies in the EPS model by u sing halo merger 
trees from iV-body simulations directly j Kauffmann et al.l 
19991; iBenson et al.ll2000l ; iHellv et al.ll2003l ; lKang et al.ll2005l ; 
Springel et ai1l2005D . As we find in this paper, obtaining ro- 
bust halo merger rates and merger trees requires rich halo 
statistics from very large cosmological simulations as well 
as careful treatments of systematic effects due to different 
algorithms used for, e.g., assigning halo masses, construct- 



ing merger trees, removing halo fragmentation events, and 
choosing time spacings between simulation outputs. 

The aim of this paper is to determine the merger rates 
of dark matter haloes as a function of halo mass, merger 
mass ratio (i.e. minor vs major), and redshift, using numeri- 
cal simulations of the ACDM cosmology. This basic quantity 
has not been thoroughly investigated until now mainly be- 
cause large catalogues of haloes from finely spaced simula- 
tion outputs are required to provide sufficient merger event 
statistics for a reliable construction of merger trees over a 
wide dynamic range in time and mass. We achieve this goal 
by using the public database of the Millennium simulation 
l|Springel et alj|2005l ), which follows the evolution of roughly 
2 x 10 7 dark matter haloes from redshift z = 127 to z = 0. 
This dataset allows us to determine the merger rates of dark 
matter haloes ranging from galaxy-mass scales of ~ 1O 12 M 
over redshifts z — to ~ 6, to cluster-mass scales up to 
~ lO 15 A/0 for z — to a few. We are also able to quantify 
the merger rates as a function of the progenitor mass ratio 
£, from major mergers (£ > 0.1) down to minor mergers of 
£ ~ 0.03 for galaxy haloes and down to £ ~ 3 x 10~ 4 for 
cluster haloes. 

The inputs needed for measuring merger rates in simu- 
lations include a catalogue of dark matter haloes and their 
masses at each redshift, and detailed information about their 
ancestry across redshifts, that is, the merger tree. Unfortu- 
nately there is not a unique way to identify haloes, assign 
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halo masses, and construct merger trees. In this paper we 
primarily consider a halo mass definition based on the stan- 
dard friends-of-friends (FOF) algorithm and briefly compare 
it with an alternative mass definition based on spherical 
overdensity. 

For the merger trees, we investigate two possible algo- 
rithms for treating events in which the particles in a given 
progenitor halo end up in more than one descendant halo 
('fragmentations'). We find that these events are common 
enough that a careful treatment is needed. In the conven- 
tional algorithm used in the literature, the progenitor halo 
is linked one-to-one to the descendant halo that has inher- 
ited the largest number of the progenitor's particles. The 
ancestry links to the other descendant haloes are severed 
(for this reason we call this scheme 'snipping'). We consider 
an alternative algorithm ('stitching') in this paper, in which 
fragmentations are assumed to be artefacts of the FOF halo 
identification scheme. We therefore choose to recombine the 
halo fragments and stitch them back into the original FOF 
halo. 

Earlier theoretical papers on merger rates either relied 
on a small sample of main haloes to estimate the overall 
redshift evolution over a limited range of halo masses, or 
were primarily concerned with the merge rs of galaxies or 
subha loes. For halo mergers, for example, iGovernato et al.l 
(QUI) studied z < 1 major mergers of galaxy-sized haloes 
in an open CDM and a tilted f2 m = 1 CDM model us- 
ing JV -body simulation s in a 100 Mpc box and 144 3 par- 
ticles. iGottlober et al.l ((2001) used a sample of ~ 4000 
haloes to study the environmental dependence of the red- 
shift evolution of th e major merger rate at z < 2 in ACDM. 
Berrier et all (|2006l ) studied major mergers of subhaloes in 
iV-body simulations in a 171 Mpc box with 512 3 parti- 
cles and the connection to the observ ed close pair count s 
of g alaxies. For g alaxy merger rates, iMurali et all (|2002l ) 
and iMaller et al.1 {2006) are based on up to ~ 500 galaxies 
formed in SPH simulat ions in ~ 50 Mpc bo xes with up to 
144 3 gas particles, while iGuo fc White! (I2007D used th e semi- 
analytical galaxy catalogue of |Pe Lucia et al.l (|2006T ) based 
on the Millennium simulation. 

This paper is organised as follows. Section [2] describes 
the dark matter haloes in the Millennium simulation ( H2.10 
and how we construct the merger trees ( i]2.20 . We then 
discuss the issue of halo fragmentation and the two methods 
('snipping' and 'stitching') used to treat these events in H2.3I 
The notation used in this paper is summarised in £12.41 

Section [3] describes how mergers are counted f i)3.1|l and 
presents four (related) statistical measures of the merger 
rate (333J- The relation between these merger rate statistics 
and the analytical merger rate based on the Extended Press- 
Schechter (EPS) model is derived in Section f3. 31 

Our main results on the merger rates computed from 
the Millennium simulation are presented in Section 2] We 
first discuss the z w results and quantify the merger rates 
as a function of the descendant halo mass and the progenitor 
mass ratios using merger trees constructed from the stitch- 
ing method ( 94.10 . The evolution of the merger rates with 
redshifts up to z ~ 6 is discussed in Section T4.2I We find a 
simple universal form for the merger rates and present an 
analytic fitting form that provides a good approximation (at 
the 10-20% level) over a wide range of parameters ( i)4.30 . 

Section [5] compares the stitching and snipping merger 



rates (3S3) and presents the key results from a number 
of tests that we have carried out to assess the robustness 
of our results. Among the tests are: time convergence and 
the dependence of the merger rates on the redshift spacing 
Az between the Millennium outputs used to construct the 
merger tree ( tj5.2[l ; how the counting of binary vs multiple 
progenitor mergers affects the merger rates f i]5.3[l ; mass non- 
conservation arising from 'diffuse' accretion in the form of 
unresolved haloes during mergers ( t|5.4p ; and how the def- 
inition of halo masses and the treatment of fragmentation 
events affect the resulting halo mass function ( (15.50 . 

In Section|Sl we discuss two theoretical frameworks that 
can be used to model halo mergers: EPS and coagulation. 
A direct comparison of our merger rates and the EPS pre- 
dictions for the Millennium ACDM model shows significant 
differences over a large range of parameter space ( i]6.1|l . Sec- 
tion [672] discusses Smoluchowski's coagulation equation and 
the connection between our merger rates and the coagula- 
tion merger kernel. 

The appendix compares a third merger tree (be- 
sides snipping and stitching) constructe d from the Millen- 
nium catalogue by the Durham group (|Bower et al.l 120061 : 
lHarker et aflbood : llleilv et alj|2003h . Two additional crite- 
ria are imposed on the subhaloes in this algorithm to reduce 
spurious linkings of FOF haloes. We find these criteria to 
result in reductions in both the major merger rates and the 
halo mass function. 

The cosmology used throughout this paper is identical 
to that used in the Millennium simulation: a ACDM model 
with Q m = 0.25, Q b = 0.045, fl A = 0.75, ft = 0.73, an 
initia l power-law index n = 1, and as = 0.9 (|Springel et al.l 
2005). Masses and lengths are quoted in units of Mq and 
Mpc without the Hubble parameter h. 



2 HALOES AND MERGER TREES IN THE 
MILLENNIUM SIMULATION 

2.1 Dark Matter Haloes 

The Millennium simulation provides the largest database to 
date for studying the merger histories of dark matter haloes 
in the ACDM cosmology. The simulation uses 2160 3 particles 
with a particle mass of 1.2 x 1O 9 M in a 685 Mpc box and 
traces the evolution of roughly 2 x 10 7 dark matte r haloes 
from redshift z = 127 to z = dSpringel et al.ll2005h . 

The haloes in the simulation are identified by group- 
ing the simulation partic les using the stan dard friends-of- 
friends algorithm (FOF: iDavis et al.1 1 19851 ) with a linking 
length of b — 0.2. Each FOF halo (henceforth referred to 
as FOF or halo) is then broken into constituent subhaloes 
by the SUBFIND algorithm, which identifies dark matter 
substructure as locally overdense regions within each FOF 
and r emoves any remainin g gravitationally unbound par- 
ticles (jSpringel et al.ll200ll ). The result is a list of disjoint 
subhaloes typically dominated by one large background host 
subhalo and a number of smaller satellite subhaloes. 

Each subhalo in the catalogue is assigned a mass given 
by the number of particles bound to the subhalo; only sub- 
haloes with more than 20 simulation particles are included 
in the database. Each FOF halo is then given two definitions 
of mass: Mfof, which counts the number of particles asso- 
ciated with the FOF group, and M200, which assumes the 
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halo is spherical and computes the virial mass within the 
radius at which the average interior density of the halo is 
200 times the mean density of the universe. Mfof includes 
background particles that are unbound by the SUBFIND al- 
gorithm so it is generally larger than the sum of the subhalo 
masses. In this paper we mainly use Mfof as it is found to 
be the more robust mass definition in our merger study. We 
discuss M200 and a number of mass conservation issues in 
Section I5T41 

2.2 Merger Tree Construction 

Merger trees of dark matter haloes in the Millennium 
database are constructed by connecting subhaloes [not the 
FOF haloes) across 64 snapshot outputs: a subhalo at a 
given output is taken to be the descendant of a progenitor 
subhalo at a prior output (i.e. higher redshift) if it contains 
the largest number of bound particles in the progenitor sub- 
halo. This procedure results in a merger tree in which each 
progenitor subhalo has a single descendant subhalo, even 
though in general, the particles in the progenitor do not 
necessarily all end up in the same descendant subhalo. 

It is worth noting that merger trees in TV-body simu- 
lations are typically constructed based on the FOF haloes 
and not on the subhaloes. The standard way of assigning 
the progenitor and descendant FOF haloes in those stud- 
ies, however, is the same as the procedure applied to the 
subhaloes in Millennium discussed above; that is, the de- 
scendant halo is the halo that inherits the most number of 
bound particles of the progenitor. As will be elaborated on 
below, we call this the 'snipping' method. 

The focus of this paper is on the merger history of the 
FOF haloes rather than the subhaloes, so we must process 
the subhalo merger tree available from the public database 
to construct a consistent merger tree for the FOF haloes. We 
consider an FOF halo A to be a descendant of an earlier FOF 
halo B if B contains a subhalo whose descendant subhalo is 
in A. Progenitor FOF haloes are said to have merged when 
all their descendant subhaloes are identified with one descen- 
dant FOF. We illustrate this process in Fig.[T]with an actual 
merger tree taken from the Millennium database. The upper 
left corner, for example, shows three FOF haloes at z = 0.24 
with masses 8.5 x 10 12 , 4 x 10 11 , and 3.8 x 10 10 M© merg- 
ing into a single FOF halo at the next Millennium output 
(z = 0.21). The largest FOF halo z = 0.24 has 7 subhaloes 
(white circles) in addition to the host (sub)halo, while each 
of the two smaller FOF haloes has only one host (sub)halo. 
For clarity, the ancestral links between subhaloes are sup- 
pressed in Fig. \T\ 

2.3 Halo Fragmentation 

Even though each subhalo in the Millennium tree, by con- 
struction, is identified with a single descendant subhalo (see 
last subsection), the resulting FOF tree can have fragmen- 
tation events in which an FOF halo is split into two (or 

1 It is common practice in the literature to call the descendant 
halo the parent halo even though the parent is formed later and, 
hence, is younger than the progenitor. We avoid this confusing 
notation throughout. 



850 




Figure 1. Example of a typical FOF merger tree extracted from 
the Millennium database. Black circles denote FOF haloes; white 
circles within black circles denote subhaloes. The radius of each 
circle is proportional to the log of the mass of the object; the 
black circles are further scaled up by a factor of 1.5 for clarity. 
(The locations of the white circles within their parent FOF haloes 
are drawn randomly.) Red circles denote fragmenting subhaloes. 
The highlighted (yellow) fragmentation event is studied in Fig. [2] 
The numbers above the haloes at z = 0.24 and to the right of the 
final descendant FOF at z = correspond to the FOF masses (in 
units of 10 10 M Q ). 
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Figure 2. Left: A closeup of the highlighted (yellow region) frag- 
mentation event in Fig. [T] The middle and right panels illustrate 
how the snipping and stitching methods handle fragmentation in 
order to assign a unique descendant halo. The blue circle (centre 
panel) shows the snipped orphan subhalo, and the yellow circle 
(right panel) shows how that subhalo is stitched. The black, white, 
and red circles are the same as in Fig. [T] 



more) descendant FOF haloes. The red circles in Fig. [T] at 
z = (0.12 : 0.09) and (0.06:0.04) illustrate two such events: 
the subhaloes of the progenitor FOF halo end up in different 
descendant FOF haloes. It is important to emphasise that 
this fragmentation issue is not unique to the use of subhaloes 
in the Millennium simulation, but rather occurs in general 
in any merger tree construction where groups of particles at 
two different redshifts must be connected. This is because 
particles in a progenitor halo rarely end up in exactly one 
descendant halo; a decision must therefore be made to select 
a unique descendant. There is not a unique way to do this, 
and we explore below two methods that we name snipping 
and stitching to handle these fragmentation events. 

Fig. [2] illustrates these two methods for the fragmenta- 
tion event shown in the highlighted (yellow) region of Fig. [1] 
The snipping method is co mmonly used in the literature 
(e.g. ISheth fc Tormen||2002l ). presumably for its simplicity. 
Fragmentation events are removed by 'snipping' the link be- 
tween the smaller descendant halo and its progenitor FOF 
halo, as shown in the middle panel of Fig. [2] The fragment- 
ing progenitor FOF halo then has only one descendant FOF 
halo. We note that this method can result in a number of 
progenitor- less orphan FOF haloes (e.g., the blue subhalo in 
Fig. El. 

In this paper we investigate an alternative method that 
we name 'stitching.' This method is motivated by our obser- 
vation that about half of the fragmented haloes in the Mil- 
lennium simulation remerge within the following 2-3 outputs 
(see below). The two fragmentation events in Fig.[T]both be- 
long to this category: the fragmented haloes at z = 0.09 and 
0.04 (red circles) are seen to have remerged by the follow- 
ing output time (z = 0.06 and 0.02). This behaviour is not 
too surprising because merging haloes oscillate in and out of 
their respective virial radii before dynamical friction brings 
them into virial e quilibrium (typically on tim escales of a 
few Gyrs; see, e.g.. iBovlan-Kolchin et al"1l2008l ). During this 
merging phase, the FOF halo finder can repeatedly disasso- 
ciate and associate the progenitor haloes, leading to spurious 
fragmentation and remerger events and inflating the merger 
rate. This behaviour needs to be taken into account before 
a robust merger rate can be obtained. 

We therefore do not count remerging fragments as 




3.58 



Figure 3. Distribution of the ratio of fragmentation to merger 
events as a function of redshift. The dotted vertical lines corre- 
spond to the redshifts of the Millennium outputs. We choose 6 
redshifts (labelled) for illustrative purposes and plot the ratio of 
the number of fragmentations to the number of mergers (filled 
circles) at each redshift. A mass ratio cutoff is applied: both the 
fragments and mergers must have mass ratios exceeding 10%. 
The line emanating from each circle then traces the evolution of 
the number of fragmentation events (the number of mergers be- 
ing held fixed), which drops as subhalo fragments remerge with 
their original FOF halo. We note that about half of the subhalo 
fragments remerge within 2-3 simulation outputs. Finally, the six 
filled circles decrease with increasing redshift, reaching ~ 40% at 
2 = but dropping to ~ 5% at high z — this is primarily due to 
the increasing Az between Millennium outputs. 



merger events in the "stitching" method. Specifically, we 
group the fragmented haloes into two categories: those that 
remerge within 3 outputs after fragmentation occurs, and 
those that do not. The fragmented haloes that remerge are 
stitched into a single FOF descendant (e.g. the yellow sub- 
halo in the right panel of Fig. [2]); those that do not remerge 
are snipped and become orphan haloes. Often the fragment 
subhaloes have become members of a new FOF group that 
is otherwise unrelated to the original FOF. In such instances 
they are removed from that group and stitched into the main 
FOF descendant A further test of the dependence of our 
results on the choice of 3 outputs is described in Section [571] 
As can be seen in Fig. the snipping method will yield 
a higher merger rate than stitching due to the remerger 
events. We quantify the relative importance of these events 
in Fig. [3] where the ratio of fragmentation events to merger 
events is seen to peak at 40% for major fragmentation events 
(defined to be fragmentations where the fragment subhalo 
carries 10% or more of the halo mass) at low- z and falls off 
at high 2 where Az is large. For the fragmentation events 
occurring at a given redshift Zf in Fig. [3] (filled circles), the 



2 There is, however, one exceptional case: if a subhalo fragment 
becomes the largest subhalo of an FOF, all subhaloes in that FOF 
are stitched into the fragment's original FOF. 



© 0000 RAS, MNRAS OOO.mfTsI 



The Nearly Universal Merger Rate 5 



drop of each curve with decreasing z tracks how many of 
them have remerged by that redshift. As noted above, we 
find that about half of the fragmented haloes remerge within 
2-3 outputs (corresponding to a fixed Azj(\ + z) as the out- 
puts are log-spaced). Given a fragmentation-to-merger ratio 
of 40%, and a remerger rate of 50%, the remerging fragments 
can impact the merger rate measurements inflating them at 
the ~20% level. 

Moreover, we find that this effect is more severe for frag- 
mentations where the mass of the fragment is small relative 
to the mass of the original parent halo (we call these mi- 
nor fragmentations) . If we consider fragmentations in which 
the subhalo fragments carry between 1% and 10% of the 
original FOF mass, the fragmentation-to-merger ratio at 
z = (z = 1.6) jumps to 57% (13.3%) vs 39% (6%) for 
major fragmentations. For very minor fragmentations (sub- 
halo fragments that carry less than 1% of the total mass) the 
fragmentation-to-merger ratios are 85% and 28% at z = 
and z — 1.6 respectively. Thus we anticipate that fragmen- 
tation events will more severely pollute the minor-merger 
regime of the merger rate statistics. 



2.4 Notation 

We apply both the stitching and snipping methods and pro- 
duce FOF merger trees from the 46 Millennium outputs that 
span z = and z = 6.2. From these trees we connect dif- 
ferent outputs and generate a catalogue of descendant FOF 
haloes at the low-z (zd) output and their associated pro- 
genitor FOF haloes at the high-z (zp) output. We refer to 
this as the zp:zn catalogue and produce a number of cata- 
logues for a variety of output spacings. The redshift spacing 
is denoted by Az — zp — zd- The Millennium outputs are 
logarithmically distributed, providing fine Az down to 0.02 
(corresponding to ~260 Myrs) near z — and larger Az at 
high redshifts, e.g., Az « 0.1 at z ~ 1 and w 0.5 at z ~ 6. 
Specifically, the lowest 10 redshift outputs are at 0.0, 0.02, 
0.04, 0.06, 0.09, 0.12, 0.14, 0.17, 0.21, and 0.24. 

For a given FOF descendant halo in a zp:zo catalogue, 
we use Mo to denote its Mfo f mass, N p to denote the num- 
ber of progenitor haloes, and Mi with i £ (1,2,..., N p ) to 
denote the rank-ordered Mfof mass of the progenitors, i.e. 
Mi ^ M2 ^ . . . Mjv p . We impose a minimum mass cutoff of 
M ^ 2 x 1O 12 M on the descendant FOF halo and a cutoff 
of Mi ^ 4.8xlO lo M on the progenitors, which corresponds 
to 40 particles and is twice the minimum halo mass in the 
Millennium database. 

For certain results reported below, we make use of three 
large mass bins: 2 x 10 12 < M < 3 x 10 13 M Q , 3 x 10 13 < 
M < 10 14 M Q , and 10 14 M Q < M , referred to as the galaxy- 
scale, group-scale, and cluster-scale bins, respectively. 



3 MERGER STATISTICS AND CONNECTION 
TO EPS 

3.1 Counting Many-to-One Mergers 

Despite the fine time spacing between Millennium's outputs, 
a non-negligible number of the descendant FOF haloes have 
more than two progenitors listed in the merger tree (i.e. 




Number of Progenitors 

Figure 4. Distribution of the number of progenitors, N p , for the 
2 = 0.06:0 merger tree. There are ~ 300, 000 descendant FOF 
haloes at redshift (black) with M )2x lO 12 Af . Of these 
~ 280, 000 have 2 x 10 12 < M < 3 X 1O 13 M (galaxy-scale; dark 
blue), ~ 16,000 have 3 X 10 13 < M Q < 1O 14 M (group scale; 
red), and ~ 5,400 have Mq ^ 1O 14 M (cluster-scale; green). 



N p > 2) . For completeness, we list in Table[T]the actual num- 
ber of merger events in the Millennium simulation available 
to us after we construct the FOF merger trees. Statistics at 
five representative redshifts are shown: z ps 0, 0.5, 1, 2, and 
3. At each z, we list separately the number of FOF haloes 
that have N p — 1,2 and > 2 progenitor haloes, for three 
separate descendant mass bins. As expected for hierarchical 
cosmological models, the halo numbers drop with increasing 
z and increasing Mo. 

Fig. [4] shows the distribution of the number of progen- 
itors, f(N p ), for the z — 0.06:0 merger tree for the same 
three mass bins. Only the stitching method is shown; the 
snipping method has a similar distribution. We find that 
(62, 22, 16)% of the haloes have N p = (1, 2, > 2) identifiable 
progenitors at z — 0.06; more than half of the FOF haloes at 
z = therefore have only one progenitor at z — 0.06 and did 
not experience a merger during this redshift interval. When 
separated into different descendant mass bins, the peak of 
f(N p ) moves to higher N p for more massive haloes. For a 
fixed {zp, zd), clusters therefore tend to have more progen- 
itors, and unlike galaxy-mass haloes, very few of the cluster 
haloes are single-progenitor events (i.e. N p = 1) 

For completeness, we include all the progenitors (above 
our minimum mass cutoff of 40 particles) in our merger rate 
statistics. Since we have no information about the order in 
which the multiple progenitors merge with one another, we 
assume that each progenitor halo Mi with i ^ 2 merges with 
Mi, the most massive progenitor, at some stage between the 
two outputs. Thus a descendant halo with N p progenitors is 
assumed to be the result of a sequence of (N p — 1) binary 
merger events, where each merger event is assigned a mass 
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N p = 1 


N p = 2 


N p > 2 


Np = 1 


N p = 2 


N p >2 


Np = 


1 N p = 2 


N p >2 


0.06:0 


188,400 


65,711 


27,939 


1,063 


2,418 


13,256 


3 


25 


5,356 


0.56:0.51 


189,351 


61,718 


22,031 


1,212 


2,468 


9,374 


6 


18 


3,014 


1.08:0.99 


145,779 


68,467 


35,426 


325 


878 


7,630 





2 


1,308 


2.07:1.91 


76,298 


52,525 


39,097 


31 


77 


2,225 








129 


3.06:2.83 


30,641 


26,675 


25,072 





i 


343 








4 



Table 1. The number of merger events in the Millennium simulation that we use to determine the merger rates. Merger trees at five 
representative redshifts are shown: z 0, 0.5, 1, 2, and 3. At each z, we list the number of FOF haloes that have a single progenitor 
halo (N p = 1, i.e., no mergers), two progenitors (N p = 2, i.e. binary mergers), and multiple progenitors (N p > 2), for three separate 
descendant mass bins: 2 X 10 12 < Mq < 3 X 1O 13 M (galaxy), 3 X 10 13 < M < 10 14 M Q (group), and Mq > 10 14 M Q (cluster). Only 
progenitor haloes with mass > 4.8 X 10 10 Mq (40 simulation particles) are counted. 



ratio 



M_ 

Mi 



(1) 



which by construction satisfies £ ^ 1. This assumption ig- 
nores the possibility that two smaller progenitor FOF haloes 
merge together before merging with the most massive pro- 
genitor. Section r5.3l describes how we have tested the validity 
of this assumption and found negligible effects as long as a 
sufficiently small Az is used. 



3.2 Definitions of Merger Rates 

In this subsection we define four related quantities that will 
be used to measure the merger rates of dark matter haloes. 
Merger rates can be measured in either per Gyr or per unit 
redshift; the two sets of quantities are related by a factor 
of dt/dz. We will present most of our results in units of per 
redshift since, as we will show below, the merger rates have 
a particularly simple form in those units. 

As a starting point, we consider the symmetric merger 

rate 



B MM > (M, M', z P :z D )dMdM' , 



(2) 



which measures the mean merger rate (i.e. the number of 
mergers per unit redshift) per unit volume between progen- 
itor FOF haloes in the mass range (M,M + dM) and (M', 
M' + dM'). We compute this quantity using merger trees 
constructed between the progenitor output redshift zp and 
the descendant output redshift zd. Note that B MM /(M, M') 
has units of [number of mergers X (Az)" 1 Mpc" 3 Mg 2 ] 
and generally depends on both zp and zp,- 

Instead of the individual progenitor masses M and M', 
it is often useful to express merger rates as a function of the 
descendant FOF mass and the mass ratio of the progenitors. 
We do this by transforming B MM > (M, M')dMdM' to 



B(M ,t,z P :z D )dM d£, 



(3) 



which measures the mean merger rate (per volume) for de- 
scendant FOF haloes in the mass range (Mo, Mo + dMo) 
at redshift zn that have progenitor FOF haloes at zp with 
mass ratio in the range of (£, , where £ = Mi jM\ , i ^ 2 

as discussed in Section r3.il The quantity B(Mo,£) therefore 
has units of [number of mergers x Az" 1 Mpc" 3 Mg 1 d£ _1 ] . 
In the mass-conserving binary limit of Mo = M + M' and 
£ = M'/M (where M' < M), B MM > and B in equations ([2j) 
and ([3]) are related by a simple transformation. In practice, 



the relation between the two quantities is complicated by 
multiple mergers and imperfect merger mass conservation. 

Since the halo abundance in a ACDM universe decreases 
with increasing halo mass, many more haloes contribute to 
the merger rates in equations ((2| and © in the lower mass 
bins of M, M', or Mo. It is useful to normalise out this effect 
and calculate the mean merger rates per halo. To do this, 
we divide out the number density of the descendant FOF 
haloes from the merger rate B and define: 



B = B(Mq,^z p :z d ) 
n ~ n(M ,z D ) 



(4) 



which measures the mean number of mergers per 
halo per unit redshift for a descendant halo of 
mass Mo with progenitor mass ratio f; the units are 
[number of mergers/number of descendants X (Az)" 1 (d£) -1 ] 
which is dimensionless. The mass function n(Mo,z)dMo 
gives the number density of the descendant FOF haloes 
with mass in the range of (Mo, Mo + dMo). 

The differential merger rates denned above can be in- 
tegrated over £ and Mo to give the mean merger rate over 
a certain range of merger mass ratios for haloes in a given 
mass range. Explicitly, the mean rate of mergers for descen- 
dant haloes in mass range Mq € [m, M] with progenitor 
mass ratios in the range £ G (x,X), 



dN n 



([m,M], [x,X],z P :z D ) 



dz 

is simply an integral over B(Mo, £, zp-.zp,) 



(•>) 



dN merge 1 

d~z = N 



B(Mo,t,z P :z D )d£dMo, (6) 



where 



nM 

N = / n(M ,z D )dM 

J m 



(7) 



is the total number of descendant haloes in the relevant mass 
range. For sufficiently small (M — m), d7V mGrgc /dz is simply 
related to the merger rate per halo, B jn, by 

/ (8, 



dz 



3.3 Connection to EPS 



The merger rates determined from the Millennium simula- 
tion can be compared to the analytic pre dictions of the Ex- 
tended Press Schech ter (EPS) formalism (|Bond et al.lll99ll ; 
lLacev fc Cold |l993) . To relate our per halo merger rate 
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B/n t o EPS, we begin with equation (2.18) of lLacev fe Cold 
(1993) for 



d 2 p 



d\nAM LC dt 



(Mi 



LC 



M^ c \t), 



(9) 



the probability that a halo of mass Mi c will merge with an- 
other halo of mass AM LC = M 2 iC - Mf c in time interval 
dt. Their notation (which we denote with superscripts 'LC') 
is related to ours by M 2 iC -»■ Mo, with M± c and AM LC 
mapped to our progenitor masses Mi and M 2 . As we will 
see below, the order is ambiguous due to an inconsistency 
in the EPS model that stems from the assumption of binary 
mergers. To relate d 2 p/dln AM LC /dt to B/n, we first mul- 
tiply it by n(Mi C ), then convert the variables to (Mo,£) 
(see below), and finally divide by n(Mo). 

Before presenting the actual equation relating the two 
rates, we note two caveats. First, in order to compute an an- 
alytical merger rate from EPS we must assume that merg- 
ers are binary and perfectly mass conserving, i.e., Mo = 
Mi + Mi in our notation. Neither assumption is strictly true 
in numerical simulations, e.g., Table Q] and Fig. H] show the 
distributions of the progenitor multiplicity N p . We defer to 
Section [5] for a detailed discussion of the tests that we have 
performed to quantify the binary nature and the degree of 
merger mass conservation in the Millennium simulation. 

Second, the EPS rate in equation <(9j is not symmetric 
in the progenitor masses M\ ° and AM LC , in contrast to 
our merger rate B MM i in equation © , which is constructed 
to be symmetric in the progenitor masses M and M'. We 
will therefore get different EPS rates depending on if Mi C 
is chosen to be the bigger or smaller progenitor. We will 
examine both options below: (A) £ = AM LC /M[ c < 1 and 
(B) £ = Mi C / AM LC sC 1. 

With these caveats in mind, we find that the per halo 
merger rate B/n corresponds to the following expression in 
the EPS model: 



B{M ,j,z) 
n(M ,z) 



2 dS c 



1 



7T dz cr(M') 



diner 



dlnM 



1 - 



2 (M ) 



-, -3/2 



<T 2 (M') 



(10) 

where M' can be the smaller progenitor, i.e., M' = M 2 = 
Mo£/(l + £) (option A), or the larger progenitor, i.e., M' = 
Mi = M /(l + £) (option B). The variable a 2 (M) is the 
variance of the linear density field smoothed with a win- 
dow function containing mass M, and 5 c (z) oc 1/D(z) is 
the standard density threshold, with D(z) being the linear 
growth factor. Note that the exponential dependence at the 
high mass end of the halo mass function has cancelled out 
on the right hand side of equation (|10|1 . Also note that both 
sides of equation (|10|) are for merger rates per redshift and 
not per time. 

We present our results for the merger rates determined 
from the Millennium simulation in the next section and com- 
pare them to the two EPS predictions in Section 16.11 



4 RESULTS 

Throughout this section, we report our results from the Mil- 
lennium merger tree where the fragmented haloes are han- 
dled with the stitching method. We find the merger rates 
given by the snipping method to agree with the stitching re- 




10 1 



10 12 10 13 10 14 
M(M ) 



1(V 



Figure 5. Symmetric merger rate B MM i of equation J2} as a 
function of progenitor masses M and M' computed from the 
z = 0.06:0 Millennium merger tree. The merger rates decrease 
from blue to red; the overlaid black lines are contours of constant 
merger rates. 



suits to within 25%. Details of the comparison are discussed 
in Section \5. II 



4.1 Merger Rates at z ~ 

Fig. [5] is a contour plot of the symmetric merger rate in 
equation B MM i (M, M' , zp : zn), calculated using the 
stitching merger tree constructed from the z = 0.06:0 Millen- 
nium outputs. Darker (bluer) regions denote higher merger 
rates, which are concentrated in the lower (M, M') corner 
because there are more low mass haloes. Minor mergers 
(off-diagonal) are more common than major mergers (along 
the diagonal). The lower left corner is blank due to our 
lower cutoff on the descendant FOF mass (~ 1000 particles; 
Mo > 2 x 10 12 Mq). The noisy nature of the upper right 
corner is due to limited merger statistics at ~ 1O 15 M0. 

As we discussed in Section 13.21 instead of progenitor 
masses M and M', it is often more illuminating to study 
merger rates as a function of the descendant FOF halo mass 
Mo and the mass ratio £ of the progenitors. This is shown in 
Fig- for the same dataset as in Fig. [5] The left panel plots 
the merger rate B(Mq, £, 0.06 : 0) of equation ([3]) against the 
progenitor mass ratio £ for fixed bins of descendant FOF 
mass Mo. We observe that the merger rate B(Mo,£) is a 
power-law in the progenitor mass ratio £ when £ < 0.1 and 
shows an upturn in the major merger regime. The power- 
law index is close to —2 and is nearly independent of the 
descendant mass Mo. More precise values are given in the 
fitting form in equation (|12p and Table [2] below. 

The main quantity we study in this paper is the mean 
merger rate per descendant halo, B/n, of equation Q, shown 
in the right panel of Fig. [S] The rising amplitude of B with 
decreasing Mo is remarkably largely removed when B/n is 
plotted: the curves in the left panel for different Mo mass 
bins collapse onto nearly a single curve in the right panel. 
This behaviour indicates that the merger rate per halo is 
nearly independent of the descendant halo mass. This weak 
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$ = M l jM l i = Mi/Mi 



Figure 6. Left panel: Mean merger rate B(Mq,£) of equation for the z = 0.06:0 merger tree as a function of the mass ratio of 
the progenitors, £, for bins of fixed descendant halo mass Mo (colour coded from black to red for increasing Mo). The overlaid dashed 
blue lines are from our fitting formula in equation 11121 . Note that the presence of a fixed minimum mass resolution (4.7 X 10 10 M Q ) 
corresponds to a minimum mass ratio £ that decreases as Mo increases. Right panel: Mean merger rate per halo, B(Mo,£)/n(Mo), of 
equation {1} for the same tree. Dividing out the halo number density n(Mo) brings the curves on the left panel to nearly a single curve, 
indicating B/n has very weak dependence on Mq. 



mass depen dence is further illust rated in Fig. [Jj and is also 
reported in iGuo fc Whitl (|2007T I. As we will quantify in 
Section [4.31 below, the dependence on Mo is approximately 
oc MS' 08 . 

Our lower cutoff of 40 particles for the progenitor FOF 
halo mass implies a lower cutoff in the mass ratio of £ ^ 4.8 x 
1O 1O M /M O . This resolution cutoff is seen in the left panel of 
Fig. [6] where we have sufficient halo statistics to measure the 
merger rates for the higher mass haloes (lower curves) down 
to very minor mergers, e.g., £ < 10 -3 for Mo > 5 x 1O 13 M ; 
whereas the dynamic range is smaller for galaxy-size haloes, 
e.g., £ > 0.01 for M < 5 x 10 12 M Q . 

The present-day merger rates shown in Fig. [5] are all 
obtained from the z = 0.06:0 merger tree. The low-redshift 
outputs available from the Millennium database in fact have 
a smaller spacing of Az ~ 0.02. We use the 0.06:0 merger 
tree to avoid any edge effects arising from our stitching crite- 
rion that only subhalo fragments that remerge within three 
outputs are stitched together (see Sec. 12. 3p . In practice, this 
precaution is not critical and we find little difference between 
the 0.06:0 and 0.02:0 results. 

4.2 Merger Rates at Higher Redshift 

Figs. HJ] and [7] summarise our results for the z = merger 
rates. At higher redshifts, the Millennium database provides 
sufficient halo statistics for us to measure merger rates up 
to 2 ~ 6. The results are shown in Fig. [5] where we plot the 
merger rate per unit time (upper panel), dN mclge /dt, and 
per unit redshift (lower panel), dN mmse /dz, as a function 
of redshift for three ranges of descendant masses (galaxy, 
group, cluster) and four ranges of progenitor mass ratios 
(£ ^ 1/3, 1/10, 1/30, and 1/100). Errors are computed as- 




2 x 10 12 10 13 10 14 10 15 

M O (M ) 

Figure 7. Mean merger rate per halo (per unit z), dN m /dz, 
as a function of descendant mass, Mo, for various ranges of the 
progenitor mass ratio £. The upper curves include increasingly 
more minor mergers. The z=0.06:0 merger tree is used. Note the 
weak mass dependence over three decades of mass. The error bars 
are computed assuming Poisson counting statistics in both the 
number of mergers and the number of haloes. 



suming Poisson statistics for the number of mergers and 
haloes. We have suppressed merger rates with poor merger 
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Lookback Time (Gyr) 
4.9 7.6 10.2 11.4 12.0 12.4 12.6 




1 + 



Figure 8. Upper panel: Mean merger rate per halo (per Gyr), dN m /dt, as a function of redshift for three bins of descendant mass Mo 
and four ranges of progenitor mass ratio £ from the Millennium simulation (using the stitching tree). The overlaid lines plot the best-fit 
power laws, (l + z)™ m , with n m ranging from 2.03 to 2.29 (labelled). Note that power laws are reasonable fits at z > 0.3 but undcrprcdict 
the Millennium rates at lower z. Lower panel: Same as the upper panel but showing the merger rate dN m /dz per unit z instead of per 
Gyr. The dotted grey lines here show our fitting formula in equation l|12|l. which is tuned to provide close fits at low z. In both panels, 
the error bars are computed assuming Poisson counting statistics in both the number of mergers and the number of haloes, and the curve 
for galaxy-scale haloes (triangles) with £ (green) is suppressed because such minor mergers fall below the simulation resolution 

limit. 



statistics (and, therefore large error bars) to keep the plots 
legible. 

The mean merger rate per Gyr (upper panel) is seen to 
increase at higher z. We have fit power laws to each Mo and 
£ range (dotted curves) of the form 



dN n 



dt 



<x (l + z) r 



(11) 



and find n m ~ 2 — 2.3 for the ranges of Mo and £ shown. 
The Millennium merger rates are seen to flatten out slightly 
at low z and deviate from a power law when the cosmolog- 



ical constant starts to dominate the energy density of the 
universe. 

A large number of merger rate statistics can be easily 
read off of Fig. [5] For example, at around z = 2 (z — 4) 
every FOF halo on average experiences ~ 2-4 (10) minor 
mergers (£ < 1/30) per Gyr, and about 10-20% (70-90%) of 
FOF haloes experience a major merger (£ > 1/3) every Gyr. 

Unlike the rising dNm CTgc /dt, the merger rate per unit 
redshift, dNm eTge /dz, shows a remarkably weak dependence 
on z in Fig. [8] (lower panel), increasing only slightly between 
z — and 1 and staying nearly constant for z > 1 for all 
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ranges of Mo and £ shown. The overlaid curves are computed 
by integrating over the fitting form for B/n to be discussed 
below (Sec. [O]). 

At z > 0, Fig. [5] shows that the dependence of 
dN meige /dz on progenitor ratio £ and descendant mass Mo 
is similar to the z — merger rates shown in Fig. [6] minor 
mergers occur more frequently than major mergers, and the 
dependence on Mo is weak, with galaxy-scale haloes (trian- 
gles) on average experiencing fewer mergers (per halo) than 
cluster-size haloes (squares). 



4.3 A Universal Fitting Form 

We now propose a fitting form that can be used to approxi- 
mate the halo merger rates in the Millennium simulation dis- 
cussed in the last two subsections to an accuracy of 10-20%. 
The key feature we will use to simplify the fit is the nearly 
universal form of the merger rate (per halo) B(Mo,£)/n 
shown in the right panel of Fig. [6] and the weak redshift 
dependence shown in the bottom panel of Fig. [8] We find 
that the following functional form works well: 



1.30 



n(M ,z) 



Mo 
M 



f exp 



d5^ v 
dz 



(12) 



where M = 1.2 x 10 12 Af Q is a constant and 8 c (z) oc \/D(z) 
is the standard density threshold normalised to S c = 1.686 
at z = 0, with D(z) being the linear growth factor. Note that 
equation (|12[) is separable with respect to the three major 
variables Mo, £, and z. 

The form of the redshift dependence in equation (|12[) 
is chosen so that r\ = 1 corresponds to the EPS predic- 
tion in equation (|10|l . In addition, this form has weak z 
dependence at z > 1 since the growth factor approaches 
that of the Einstein-de Sitter model, S c (z) = 1.68(1 + z), 
and dS c /dz approaches a constant. This behaviour matches 
the weak redshift dependence seen in the Millennium merger 
rate (bottom panel of Fig. 

To determine the parameters in equation (|12p . we fit 
simultaneously to all redshifts z < 1, mass ratios £ > 10~ 3 , 
and masses 10 12 < M < 10 14 M Q . The B/n data points 
are weighted using their Poisson distributed errors. The re- 
sulting fits are plotted as dotted curves in Figs. [6] and [8] 
and the fitting parameters are given in Tabled along with 
the overall reduced xi obtained by fitting to all redshifts 
z < 1 simultaneously. In addition to computing a global xt 
we also compute a local xl ( z ) a t each redshift and find rel- 
atively good convergence across the z < 1 redshift range: 
xt(z) remains below 1.5 for stitching and below 2 to 3 for 
snipping. 

We note that the fitting form of equation ()12|) does not 
appear symmetric in the progenitor masses Mi and M2 be- 
cause by construction, £ = M2/M1 < 1. However, for any 
pair of progenitors, we identify Mi with the more massive 
and M2 with the less massive progenitor and then compute 
£ = M2/M1 < 1. This procedure yields the same £ and 
therefore the same B/n regardless of the order of the in- 
put progenitors, in contrast to the EPS model discussed in 
Section [3731 which is intrinsically asymmetric in Mi and M2. 




Mi/Mi 



Figure 9. The ratio of the snipping and stitching B/n as a func- 
tion of mass ratio £ computed using the 0.06:0 catalogue for a 
variety of mass bins in the range 2.4 X 1O 12 M0 (black) $C M ^, 
1.3 X 10 14 (red). We find differences at the 25% level at low £ with 
the snipping method consistently predicting a higher merger rate 
at all £. We attribute this to the population of remerging orphan 
haloes. 



5 TESTS 

5.1 Snipping vs Stitching Trees 

Fig. [9] shows the ratio of the z = per-halo merger rate 
B/n from the snipping and stitching methods. Overall, the 
merger rates given by the two methods differ by no more 
than 25% over 2-3 orders of magnitude in both the pro- 
genitor mass ratio £ and the descendant mass Mo. Within 
this difference, however, Fig. [9] and Table [2] show that the 
snipping method systematically yields a higher merger rate 
and a steeper slope in the ^-dependence than the stitching 
method. These additional merger events come from the or- 
phaned subhaloes that are first snipped and subsequently 
remerge (see Fig. Moreover, as discussed in Section [2.31 
the fragmentation-to-merger ratio is higher for more minor 
subhalo fragments (those with low fragment-to-FOF mass 
ratios). There are therefore more remerging orphan haloes 
with lower £, leading to the larger difference between snip- 
ping and stitching at low £ seen in Fig. [9] 

A remaining issue is our choice of the stitching criterion: 
as described in Section 12.31 we stitch only FOF fragments 
that are observed to remerge within the next 3 outputs. This 
choice is motivated by the fact that about half of the halo 
fragments at a given output will have remerged within three 
outputs (see Fig. [3)l , and that such a small Az criterion will 
allow us to effectively compute instantaneous merger rates. 
We have tested this criterion further by implementing a more 
aggressive stitching algorithm that stitches all fragments, 
regardless of whether they eventually remerge. We call this 
oo-stitching. This algorithm represents the opposite limit 
to the snipping method and may err on the side of under- 
estimating the merger rates since it would stitch together 
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Method A g a (3 7 n \t 

Snip 0.0101 0.017 0.089 -2.17 0.316 0.325 1.86 
Stitch 0.0289 0.098 0.083 -2.01 0.409 0.371 1.05 



Table 2. Best fit parameters for equation (1 1 2 D . 



close-encounter fly-by events that do not result in actual 
mergers within a Hubble time. We find the amplitude of B /n 
from oo-stitching to be lower than that from the 3-stitching 
by up to ~ 25%, similar in magnitude but opposite in sign 
to the difference between snipping vs 3-stitching shown in 
Fig. [9] The fitting form in equation fjl2p works well for oo- 
stitching, where the best fit parameters are A = 0.0344, £ = 
0.125, a = 0.118,/? = -1.921,7 = 0.399, and n = 0.853. 
This algorithm shows excellent convergence properties (see 
i|5.2p and excellent mass conservation properties f i]5.4[l but 
alters the FOF mass function by a few percent. 

Since the snipping algorithm tends to inflate the merger 
rate and the oo-stitching algorithm tends to under-estimate 
it, we believe the 3-stitching used in all the results in Sec- 
tion [4] should be a fairly robust scheme. 

5.2 Convergence With Respect to Az 

We have performed a number of tests to quantify the depen- 
dence of our merger rate results on the choice of Az between 
the Millennium outputs used to construct the merger trees. 
It is not a priori clear which value of Az is optimal: small Az 
can result in poor merger statistics since most haloes would 
not have had time to merge; whereas large Az does not have 
the time resolution to track individual merger events accu- 
rately and also runs the risk of smearing out real redshift 
dependent effects. The optimal Az may also vary with red- 
shift. 

Our first test focuses on z ~ mergers and quantifies 
how B/n varies with the Az used to construct the trees. 
Fig. 1101 shows the ratios of B/n for five pairs of progeni- 
tor and descendant redshifts: (z P ,z D ) =(0.02:0), (0.04:0), 
(0.06:0), (0.12:0), and (0.24:0), corresponding to a time in- 
terval of At = 0.26,0.54,0.83,1.44, and 2.77 Gyr, respec- 
tively. For the stitching method (left panel), there is excel- 
lent convergence for Az < 0.12 (panels A-F), where the ra- 
tios of B/n are centred around 1 and rarely deviate beyond 
the 10% level (dotted line). For Az = 0.24 (panels G-J), the 
ratios start to drop below unity. This is consistent with the 
slowly rising merger rates with increasing z shown in Fig. [8] 
Thus, the stitching method yields merger trees with robust 
Az convergence properties near z — 0, and we have chosen 
Az — 0.06 to compute the merger rates in earlier sections. 

The snipping method (right panel) shows inferior Az 
convergence. The B/n computed with smaller Az consis- 
tently show higher merger rates than those computed with 
larger Az. Moving up the left column (panels G,D,B,A), 
we observe only some degree of convergence. Better conver- 
gence is seen along the main diagonal (panels A,C,F,J) in 
order of increasing Az. In particular, panels C and F show 
excellent convergence properties (to the 10% level) centred 
around (0.06:0). To emphasise that the problem is with Az 
and not with a particular output (say, any possible edge 
effects at z = or 0.02), we show in panels A',B',C the 
ratios of B/n computed using three merger trees with the 



same Az — 0.02 but centred at progressively higher z: z = 
(0.02:0), (0.04:0.02), and (0.06:0.04). The agreement is ex- 
cellent, in striking contrast to panel B. Based on these tests, 
we have chosen to use Az = 0.06 for the snipping method. 

We believe that the snipping method has inferior Az 
convergence properties because of the remerging orphan sub- 
haloes (see Section 12.31 and Fig. [2}. These fragmentation 
events are sewn together in the stitching scheme and there- 
fore do not contribute to the merger rates. In the snipping 
scheme, however, the snipped events provide a fresh supply 
of haloes, many of which remerge in the next few outputs. 
This effect artificially boosts the merger rate across small 
Az. 

Our second Az convergence test is performed at all red- 
shifts. We test three types of spacing: (1) Adjacent spac- 
ing uses adjacent catalogues, e.g., at low z, it uses (0.02:0), 
(0.04:0.02), (0.06:0.04); (2) Skip 1 spacing skips an output, 
e.g., (0.04:0), (0.06:0.02)...., and (3) Skip 2 spacing skips two 
outputs, e.g., (0.06:0), (0.09:0.02), and so on. Fig. ITT1 shows 
dNmerge/dz computed using these three Az for galaxy-mass 
haloes. We again see excellent Az-convergence for the stitch- 
ing method at z < 1.5 (left panel) and worse Az-convergence 
for the snipping method (right panel) . The latter follows the 
behaviour seen in Fig 1101 with adjacent spacing (Az = 0.02 
at z = 0) over-predicting the merger rate. 

At higher redshifts (z > 1.5), Fig. 1111 shows that the 
merger rates in the minor merger regime differ by up to 
~ 15% depending on which of the three types of spacing is 
used. This difference is not likely to be due to the fragmenta- 
tion events since as Fig. 0shows, the ratio of fragmentation 
to merger events is 40% near z = but drops to <10% for 
z > 3. Rather, we believe that the inferior Az convergence 
at high z is due to the increasing Az between Millennium 
outputs (e.g., the smallest Az is ~ 0.5 at z ~ 6 vs Az = 0.02 
at z « 0) and the inaccuracy of the multiple counting order- 
ing assumption for large A2: (see Section f5 . 3 f) . At high z, we 
therefore advocate using the finest output spacings available 
in the Millennium database, noting the good time conver- 
gence for major mergers (£ > 1/3) but ~ 15% variations in 
the minor merger rates. 



5.3 Multiple vs Binary Counting 

As discussed in Section T3. II for descendant FOF haloes with 
more than two progenitor haloes, we include all progeni- 
tors when we calculate the merger rates for completeness. 
Since mergers are often assumed to be binary events, we have 
tested the difference between our multiple counting results 
and those obtained by counting only the two most massive 
progenitors of a given descendant halo. Fig.[T2]compares the 
merger rates (per halo), B(Mo,£)/n, for these two counting 
methods (dashed: multiple; solid: binary) as a function of 
the progenitor mass ratio £ for three descendant masses Mq 
(increasing from left to right). In each panel, results from 
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Figure 10. Az Convergence Matrix (stitching left, snipping right). Each subplot is the ratio of B/n for two different catalogues (labelled). 
The dashed lines denote equality and the dotted lines are the 10% deviation levels. The ratios are presented for a variety of mass bins 
with the high mass bin highlighted in thick blue (or red) and the low mass bin highlighted in thick black. 




Figure 12. Comparisons of merger rate per halo, B/n, computed via multiple counting (dashed lines) and binary counting (solid lines) 
for 4 merger trees with increasing Az. Three descendant mass bins are shown (from left to right). We find the multiple counting rate 
to be in excellent agreement regardless of Az of the tree, whereas the binary counting B/n curves fall off from the observed power-law 
behaviour towards lower £. 



four merger trees using increasing Az of 0.02, 0.06, 0.12, 
and 0.24 are shown. 

The most notable trend in Fig. [T^] is that the multiple 
counting method gives similar merger rates regardless of Az, 
indicating good time convergence in the results (as we have 
discussed in detail in Section l5T2"jl . The binary counting rates, 
on the other hand, deviate increasingly from the multiple 
rates when larger Az are used because the binary assump- 
tion becomes less valid for larger Az. As a function of £, 
the binary and multiple merger rates match well in the ma- 
jor merger regime but deviate significantly for small £. This 
occurs because binary counting counts only the two most 
massive progenitors and ignores the additional (typically 
low-mass) progenitors. It therefore closely approximates the 



major-merger rates but under-estimates the minor-merger 
regime of the multiple counting result. 

Fig. [12] suggests that for a given minimum mass reso- 
lution (i.e. a minimum £), there is a corresponding Az for 
which the binary counting method is a good approximation. 
For example, for 6 x 1O 13 M0 haloes (centre panel), the bi- 
nary and multiple merger rates are similar down to £ ~ 0.05 
for Az ~ 0.12, and down to £ ~ 0.005 when Az is decreased 
to 0.02. Thus the multiple counting B/n can be thought of 
as the small- Az limit of the binary B/n. 

Another test we have performed is on the ordering of 
mergers assumed in the multiple counting method described 
in Section T3. II There, we assumed that the less massive pro- 
genitors M2, M3, ... each merged with the most massive pro- 
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Figure 11. Merger rate d/V mcrgG /eiz computed using three types 
of redshift spacings: adjacent, skip 1 and skip 2 (see text). For 
clarity, only the galaxy-mass haloes are shown; the group and 
cluster haloes behave similarly. 



genitor Mi and not with one another. This assumption is 
motivated by the fact that satellite haloes in iV-body cos- 
mological simulations are typically seen to accrete onto a 
much more massive host halo as a minor merger event in- 
stead of merging with another satellite halo. We have quan- 
tified the validity of this assumption by taking large Az in 
the Millennium outputs, applying this ordering, and check- 
ing against the actual merging order among the progenitors 
when finer Az is used. (Of course, we cannot do this for the 
minimum Az = 0.02 available in the database.) The frac- 
tion of misordering naturally rises with increasing Az due 
to the degraded time resolution, but for Az < 0.06, we find 
the fraction of progenitors to have merged with a progenitor 
other than Mi to be < 10%. Most of the mergers among 
multiple progenitors, therefore, do occur between the most 
massive progenitor and a less massive progenitor, as we have 
assumed. 



5.4 Mass Conservation and "Diffuse" Accretion 

Thus far we have analysed mergers in terms of the progenitor 
halo mass M; and the descendant halo mass Mq. Mergers 
are, however, messy events, and the sum of Mi does not 
necessarily equal Mo. To quantify this effect, we define a 
"diffuse" component, AM, for a given descendant halo: 



Mo ]>\\/. • AM . 



(13) 



where AM is diffuse in the sense that it is not resolved as 
distinct haloes in the simulation. A non-zero AM can be 
due to physical processes such as tidal stripping and diffuse 
mass accretion that cause a net loss or gain in halo mass 
after a merger event. In simulations, additional numerical 
factors also contribute to AM due to different algorithms 
used in, for example, defining halo mass (FOF vs spherical 



AM/M 





AM/Mq 



Figure 13. Distributions of AM/Mq from the z = 0.06:0 Mil- 
lennium merger tree (using stitching; snipping is nearly identical) 
computed using the Mfof mass (top panel) vs A/200 virial mass 
(bottom panel). The solid vertical line is the median of the distri- 
bution for the galaxy-mass bin and the dashed line is the mean. 
We note a longer negative AM tail for the M200 tree when com- 
pared to the MpQF tree. Note, however, that the peaks of the 
two distributions are in good agreement (AM/Mq ~ 2.5%) 



overdensity) . AM therefore does not necessarily have to be 
positive in every merger event. 

Fig. [T3] shows the distribution of AM /Mo for the 
z =(0.06:0) Millennium merger tree. Only haloes that have 
experienced mergers between these two redshifts (i.e. those 
with more than one progenitor) are plotted. The snipping 
tree (not shown) gives a very similar distribution as the 
stitching tree shown here. For comparison, we have com- 
puted AM/Mq using the two different halo mass definitions 
Mfof and M2oo- The distribution shows a prominent peak 
at AM/Mq ~ 2.5 % for both Mfof and M 20 o, indicating 
that in the majority of the merger events between z = 0.06 
and 0.0, ~ 97.5% of the mass of the descendant halo comes 
from resolvable progenitor haloes. 

Even though the two distributions in Fig. 1 131 have simi- 
lar peaks, the M200 mass definition produces longer AM/Mq 
tails than the Mfof mass definition, and the mean of the 
distribution (dotted vertical line) is negative for M2oo- We 
believe this is because mass definitions based on the assump- 
tion of spherical symmetry (as M200 does) have difficulties 
assigning accurate mass to non-spherical FOF haloes and 
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Figure 14. Ratios of the Millennium halo mass function (com- 
puted from the stitching trees) to the fit of ljenkins et al.l (2001) 
using the Mpo F mass. The results for the snipping trees are vir- 
tually identical. We note a significant deviation of ~ 25% at z ~ 3 
between the Jenkins fit and the Millennium mass function. 

tend to underestimate the halo mass (see, e.g.. lWhitdl200lh . 
Mfof, on the other hand, can account for all the mass in 
a given FOF object that is identified as 'merged' by the 
FOF halo finder well before virialization. As discussed in 
Section \2. II we have been using the Mfof mass thus far. 

Our main results on merger rates in Section [4] were de- 
termined for numerically resolved dark matter haloes; they 
are therefore not affected directly by the fact that AM is 
generally non-zero for merger events. We find, however, that 
AAf/A/o increases with Az, and this diffuse accretion com- 
ponent makes an important contribution to the mass growth 
of a halo over its lifetime. We will explore the growth of 
haloes in further detail in subsequent papers. 

5.5 Halo Mass Function 

The mass function of dark matter haloes in principle de- 
pends on both the definition of halo mass and the algorithm 
used to treat fragmentation events. As illustrated in Fig. [5] 
the snipping method by construction preserves the origi- 
nal FOF mass function, while the stitching scheme modifies 
it slightly as it rearranges fragmented subhaloes between 
FOFs. We find that the impact on the mass function is neg- 
ligible (less than ~ 0.25%) so will use the stitching result 
below. 

Fig. [14] shows the ratio of the Millenn ium halo mass 
function to the fits by I Jenkins et all (|200ll ) for Mfof at 
four redshifts: z ~ 0, 1, 2, and 3. The fit of Jenkins et al 
is accurate to better than 10% for low redshift (z < 1), 
but it underestimates the Millennium halo abundance by 
> 25% at the high mass end for z > 1. This discrepancy 
is present but not o b vious on the log-lo g plot in Fig. 2 of 
ISpringel et"al1 ([20051 ). iLukicet al.l i[2007l ) also noted this dif- 
ference. Since the stitching and snipping mass functions are 
virtually identical, this appears to be a discrepancy between 



the Millennium FOF catalogue and the fit of ljenkins et al.l 
(|200ll ). 



6 THEORETICAL MODELS FOR MERGER 
RATES 

6.1 Extended Press-Schechter Model 

In Section 13.31 we discussed how our merger rates are re- 
lated to the conditional probabilities in the EPS model and 
obtained equation (|10p . where there are two choices for the 
definition of the progenitor mass M' since the EPS model 
is not symmetric in the two progenitor masses. In Fig. [15] 
we show the ratio of the EPS prediction to our Millennium 
B/n at z — 0, where we have computed the EPS rates given 
by the right-hand-side of equation (|10[) using the same cos- 
mological parameters as for the Millennium simulation. We 
compute the variance of the smoothed linear density field, 
a 2 (M), in the ACDM cosmology using the power spectrum 
fit provided in lEisenstein fc Hul(ll999l ). 

Fig. ll5l shows that EPS underpredicts the z — rate for 
minor mergers by up to a factor of ~ 5, and overpredicts the 
rate at £ > 0.05, indicating that the dependence of the EPS 
merger rate on £ is shallower than our B/n ~ , where the 
best-fit (5 is —2.17 and —2.01 for the snipping and stitching 
methods, respectively (see Table [2j. In terms of the descen- 
dant mass Mo, the dependence of the EPS rate is too steep 
compared to our B/n, leading to the spread in each bun- 
dle of curves in Fig. 1151 The two choices of M' in EPS are 
seen to lead to different predictions. Assigning M' to be the 
smaller progenitor (option A) results in a somewhat smaller 
discrepancy than option B. 

Fig- HS]compares the rates at z = 0. At higher redshifts, 
we find the Millennium merger rate to evolve as oc (dS c /dz) v , 
where n ~ 0.37 (see equation [12] and Table [2} and is shal- 
lower than the EPS prediction of r\ — 1 in equation (|10p . 
Since the functional forms of both our fit for B/n and the 
EPS expression are separable with respect to Mo, £ and z, 
the z = curves in Fig. [TJlwill maintain the same shape at 
higher z, but the amplitude of the ratio will increase. For 
instance, the ratio shown in Fig. [15] will be increased by a 
factor of 1.26, 1.32, 1.34, and 1.35 at z = 1, 2, 4, and 6 re- 
spectively. The discrepancy between the Millennium results 
and the EPS predictions is therefore even worse at higher z. 

Given that the Press-Schechter mass function is known 
not to match the halo abundances in simulations very 
closely, it is not particularly surprising that the EPS merger 
rates in Fig.[T3]do not match the Millennium results closely. 
The substantial discrepancy, however, does highlight the 
limitation of the EPS model and provides the motivation 
to build more accurate merger rates based on improved 
PS mass functions. We address this issue in separate pa- 
pers l|Zhang et al.l [2008) . in which we investigate a moving 
density-barrier algorithm to generate merger trees that pro- 
duces a better match to simulation results than the constant 
barrier of the PS model. 

6.2 Halo Coagulation 

The merging of dark matter haloes is, in principle, a coagu- 
lation process. Coagulation is often modelled by the Smolu- 
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Figure 15. Comparison between our Millennium merger rate 
(from the fit) and the two predictions of the Extended Press- 
Schechter model. The ratio of B/n from EPS to Millennium is 
plotted. Blue and red label the two options in assigning progen- 
itor masses in the EPS model (see text); within each colour, the 
set of curves from bottom to top denotes increasing Mo bins, 
from ~ 10 12 A/ Q to ~ 3 X 1O 14 M . The EPS model is seen to 
overpredict the major merger rate by up to a factor of ~ 2 and 
undcrprcdicts the minor merger rate by up to a factor of ~ 5. 



chowski coagulation equation l|Smoluchowskil [l916h . which 
governs the time evolution of the mass function n(M, t) of 
the objects of interest with a coagulation kernel. In the ab- 
sence of fragmentations, the time change of n is given by 



dn(M) 
di 



A(M', M-M')n(M')n(M-M')dM' 



A(M, M')n(M')n(M)dM' , 



(14) 



where the first term on the right-hand side is a source term 
due to mergers of two smaller haloes of mass M' and M — M', 
while the second term is a sink term due to haloes in the 
mass bin of interest merging with another halo of mass M', 
forming a halo of higher mass M+M' . When applied to hier- 
archical structure formation, A(M,M'), the symmetric co- 
agulation kernel (in units of volume/time), tracks the prob- 
ability for a halo of mass M to merge with a halo of mass 
M . Our merger rate per halo, B/n, can be simply related 
to A by 



A(M, M') 



B(M,M') 
n(M)n(M') 



(15) 



We note, however, that the coagulation equation in the form 
of equation (|14[) is valid only for mass-conserving binary 
mergers. As seen throughout this paper, these assumptions 
are not strictly true in numerical simulations, and modi- 
fications are required to account for the issues that have 
been discussed thus far, such as net mass gain or loss (i.e. 
AM ^ 0), multiple merger events, and halo fragmentation. 
While the relative errors may be small when integrated over 
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Figure 16. Same as Fig. 1151 only now comparing the merger 
rates from Benson et al. (2005) to the Millennium B/n for initial 
power spectrum index n = —2 (blue) and n = —1 (red). 



a small time interval, repeated application of equation (|14[) 
using equation (|15[) may not yield robust results. 

Assuming that n(M) is the Press-Schechter mass func- 
tion, iBenson et all (2005) have developed numerical tech- 
niques to construct the coagulation kernel for self-similar 
cosmological models with initial power-law power spectrum 
P(k) oc k n . Their technique is underconstrained and does 
not yield a unique expression for A(M, M'). In order to pick 
out a particular solution, a regularisation condition was ap- 
plied to force A{M, M') to vary smoothly. We have trans- 
formed the coordinates of their fits to A(M, M') to compare 
their results with our merger rate B/n. Fig. [16] shows the ra- 
tio of their fitting formula to the Millennium ACDM merger 
rate for spectral indices n = — 1 and —2 as a function of pro- 
genitor mass ratio £ for various descendant halo mass bins. 
The difference can be up to a factor of several. 



7 CONCLUSIONS AND DISCUSSIONS 

In this paper we have computed the merger rates of dark 
matter FOF haloes as a function of descendant halo mass 
Mo, progenitor mass ratio f, and redshift z using the merger 
trees that we constructed from the halo catalogue of the 
Millennium simulation. Our main results are presented in 
Figs. [6] to [8] which show very simple and nearly separable 
dependence on Mo, £, and z. The mean merger rate per 
descendant FOF halo, B/n, is seen to depend very weakly 
on the halo mass Mo (Fig. [6] right panel and Fig. [TJl . As 
a function of redshift z, the per halo merger rate in units 
of per Gyr increases as (1 + z) a , where a ~ 2 to 2.3 (top 
panel of Fig. |HJ, but when expressed in units of per red- 
shift, the merger rate depends very weakly on z (bottom 
panel of Fig.[SJ). Regardless of Mo and z, the dependence of 
B/n on the progenitor mass ratio, £ = Mi /Mi, is a power 
law to a good approximation in the minor merger regime 
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(£ < 0.1) and shows an upturn in the major merger regime 
(Fig. [6]l . These simple behaviours have allowed us to propose 
a universal fitting formula in equation (|12|) that is valid for 



10 12 < M < 10 15 M Q , £ > lO- 3 , and up to z ~ 6. 

Throughout the paper we have emphasised and quanti- 
fied the effects on the merger rates due to events in which a 
progenitor halo fragments into multiple descendant haloes. 
We have shown that the method commonly used to re- 
move these fragmented haloes in merger trees - the snip- 
ping method - has relatively poor Az-convergence (Figs. [101 
and Hip . Our alternative approach - the stitching method 
- performs well with regards to this issue without dras- 
tically modifying the mass conservation properties or the 
mass function of the Millennium FOF catalogue (Figs. 1131 
and []g| . 

We have computed the two predictions for merger rates 
from the analytical EPS model for the same ACDM model 
used in the Millennium simulation. At z = 0, we find the 
EPS major merger rates to be too high by 50-100% (de- 
pending on halo mass) and the minor merger rates to be 
too low by up to a factor of 2-5 (Fig. 1151) . The discrepancy 
increases at higher z. 

The coagulation equation offers an alternative theoret- 
ical framework for modelling the mergers of dark matter 
haloes. We have discussed how our merger rate is related to 
the coagulation merger kernel in theory. In practice, how- 
ever, we find that mergers in simulations are not always 
mass-conserving binary events, as assumed in the standard 
coagulation form given by equation (|14[l . Equation (|14[) will 
therefore have to be modified before it can be used to model 
mer gers in simulations. 

iGottlober et al.l (|200ll ) studied the rate of major merg- 
ers (defined to be £ ^ 1/3 in our notation) in A-body simula- 
tions and found a steeper power law dependence of oc (1+z) 3 
(at z < 2) for the merger rate per Gyr than ours. Their simu- 
lations did not have sufficient mass resolution to determine 
the rate at z > 2. It is important to note, however, that 
our B/n at redshift z measures the instantaneous rate of 
mergers during a small Az interval at that redshift. By con- 
trast, they studied the merging history of present-day haloes 
and measured only the rate of major mergers for the most 
massive progenitor at redshift z of a z = halo (see their 
paragraph 4, section 2). A detailed comparison is outside 
the interest of this paper. 

Mergers of dark matter haloes are related to but not 
identical to mergers of galaxies. It typically takes the stellar 
component of an infalling galaxy extra time to merge with 
a central galaxy in a group or cluster after their respective 
dark matter haloes have been tagged as merged by the FOF 
algorithm. This time delay is governed by the dynamical 
friction timescale for the galaxies to lose orbital energy and 
momentum, and it depends on the mass ratios of the gala xies 
and the orbital parameters dBovlan-Kolchin et al ] |2008l and 
references therein). In addition to this difference in merger 
timescale, the growth in the stellar mass of a galaxy is not 
always proportional to the growth in its dark matter halo 
mass. A recent an alysis of the galaxy ca talogue in the Millen- 
nium simulation (iGuo fc White! 120071 ) finds galaxy growth 
via major mergers to depend strongly on stellar mass, where 
mergers are more important in the buildup of stellar masses 
in massive galaxies while star formation is more important in 
galaxies smaller than the Milky Way. Extending the analysis 



of this paper to the mergers of subhaloes in the Millennium 
simulation will provide the essential link between their and 
our results. 

For similar reasons, our results for the evolution of the 
dark matter halo merger rate per Gyr ((1 + z) Um with 
n m ~ 2 — 2.3) cannot be trivially connected to the ob- 
served merger rate of galaxies. It is nonetheless interesting 
to note that a broad disagreement persists in the observa- 
tional literature of galaxy merger rates. The reported power 
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in N-body simulations and provided a m ore detaile d com - 
parison with recent observations by, e.g., iLin et all (|2004l ) 
that find n m < 1. They attributed such a weak redshift 
evolution in the number of close companions per galaxy to 
the fact that the high merger rate per halo at early times is 
counteracted by a decrease in the number of haloes massive 
enough to host a galaxy pair. 

The merger rates in this paper are global averages over 
all halo environments. The rich statistics in the Millennium 
simulation allow for an in-depth analysis of the environmen- 
tal dependence of dark matter halo merger rates, which we 
will report in a subsequent paper (Fakhouri & Ma, in prepa- 
ration) . 
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APPENDIX: THE DURHAM TREE 

In this paper we have used two methods to handle fragmen- 
tation events in the Millennium FOF merger trees: snipping 
and stitching. Here we disc uss and compare a third method 
used by the Durham gr oup (|Bower et al.ll2006l ; iHarker et al.l 
120061 ; iHellv et al J 120031 ) . 

The Durham algorithm is designed to reduce spurious 
linkings of FOF haloes in low-density regions. Before con- 
structing the FOF merger tree, they filter through the Mil- 
lennium FOF and subhalo database, and split up a subhalo 
from its FOF halo if (1) the subhalo's centre is outside twice 
the half mass radius of the FOF halo, or (2) the subhalo has 
retained more than 75% of the mass it had at the last out- 
put t ime at which it was an independent halo (|Harker et al.l 
2006). Condition (1) is effectively a spatial cut, while (2) is 
based on the argument that less massive subhaloes are ex- 
pected to undergo significant stripping as they merge with 
more massive haloes. This algorithm then discards the sub- 
haloes that are split off from FOF groups at z = 0, along 
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Figure 17. The ratio of the Durham merger rate B/n to our 
stitching rate B/n (section 14. 3j l as a function of progenitor mass 
ratio £ for a number of descendant mass bins ranging from ~ 
2 X 10 12 M Q (black) to ~ 1O 14 M (red). The Durham merger 
rate tends to be lower than the stitching merger rate, and suffers 
a sudden drop in the major merger regime (£ > 0.3). 



with any associated progenitor subhaloes. Around 15% of 
the original FOF haloes are split in this algorithm. 

The Durham algorithm tends to reduce the number of 
fragmented haloes in the resulting trees, but it does not 
eliminate all such events. A method much like our snipping 
method is used to treat the remaining fragmentation events. 
The resulting FOF tree is available at the Millennium public 
database along with the original Millennium tree. 

To compare with our stitching and snipping trees, we 
have repeated all of our merger rate calculations and tests 
using the Durham tree. Fig. H7l shows the ratio of the result- 
ing Durham merger rate, B/n, to that from our stitching 
tree at z = 0. The Durham rate is generally lower than our 
rate for minor mergers (by up to ~ 30%), and it drops pre- 
cipitously for major mergers (£ > 0.3). The two additional 
conditions applied to split up subhaloes in the Durham algo- 
rithm therefore appear to have eliminated most of the major 
merger events. 

Moreover, these splitting conditions in the Durham al- 
gorithm also modify the halo mass function. Fig. 1181 shows 
the ratio of the Durham mass function to the fit of Jenkins 
et al. (2001) at z = 0, 0.5, and 1 (thick solid curves with 
error bars). The ratio of our stitching mass function to the 
same fit is overlaid for comparison (thin dotted curves) . The 
Durham mass function is systematically lower: the number 
of z = haloes with M > 10 14 M Q is ~ 25% lower, and 
the difference increases at z ~ 1, affecting the halo mass 
function even at M ~ 2 x 1O 12 M . 

We believe that the deficit of major merger events and 
massive haloes in the Durham catalogue is partially due to 
their second criterion that splits off subhaloes that have re- 
tained 75% of their original mass. This condition may indeed 
remove spurious FOF linkings in the minor merger regime, 
but major merger events tend to preserve much of the origi- 
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Figure 18. The ratio of the Durham halo mass function to the 
fit of Jenkins et al. (2001) at redshifts 0, 0.5, and 1 (thick solid 
curves with error bars) . The ratio of the halo mass function from 
our stitching method to the same fit is shown for comparison 
(thin dashed curves). The Durham algorithm tends to reduce the 
masses of massive haloes, leading to a deficit that grows to ~ 50% 
at ~ 1O 15 M0 and at higher z. 



nal progenitor masses and have been systematically split by 
the Durham algorithm. 

Finally, Fig. \T§\ (right panel) shows that the Durham 
tree has similar mass conservation properties as our stitch- 
ing tree in Fig. 1131 The distribution of the mass in the 'dif- 
fuse' component, AM /Mo, has a very similar peak of ~ 3%, 
although the negative AM events have been suppressed. For 
Az convergence (left panels; cf section [5T2|| . the Durham tree 
performs well in the minor merger regime but is consistently 
poor for major mergers, again probably due to the splitting 
condition (2) above. 

The Durham algorithm is tuned to address questions 
of galaxy evolution. The issues we have uncovered regarding 
this algorithm are specifically for the mergers of dark matter 
haloes, the subject of this paper; issues with the mergers of 
galaxies will require a separate study. 
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Note the poor convergence properties of the major merger end 
K ;^0.3). This corresponds to the region of the largest difference 
between the stitching and Durham merger rates (see Fig. I17H . 
Right panel: The distribution of AM /Mo for the z = 0.06:0 
Durham catalogue (similar to Fig. 1131 ). 
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